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The heat transport in a system of 5=1/2 large-J Heisenberg spin chains, describing closely 
Sr 2 Cu 03 and SrCu02 cuprates, is studied theoretically at T ^ J by considering interactions of 
the bosonized spin excitations with optical phonons and defects. Treating rigorously the multi¬ 
boson processes, we derive a microscopic spin-phonon scattering rate that adheres to an intuitive 
picture of phonons acting as thermally populated defects for the fast spin excitations. The mean- 
free path of the latter exhibits a distinctive T-dependence reflecting a critical nature of spin chains 
and gives a close description of experiments. By the naturalness criterion of realistically small spin- 
phonon interaction, our approach stands out from previous considerations that require large coupling 
constants to explain the data and thus imply a spin-Peierls transition, absent in real materials. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.40.Gb, 72.20.Pa 


The one-dimensional (ID) spin chains are among the 
first strongly-interacting quantum many-body systems 
ever studied mm and they remain a fertile ground for 
new ideas 0 and for developments of advanced theoret¬ 
ical and numerical mm methods. A number of physical 
realizations of spin-chain materials [MIQ] have allowed 
for an unprecedentedly comprehensive comparisons be¬ 
tween theory, numerical approaches, and experimental 
data [IIHS]. Current theoretical challenges for these 
systems include their dynamical, non-equilibrium, and 
transport properties |14H20] . The transport phenomena 
are particularly challenging as the couplings to phonons 
and impurities, perturbations that are extrinsic to the 
often integrable spin systems, become crucial 

In this Letter, we address the problem of heat trans¬ 
port in ID spin-chain systems by considering coupling 
of spins to optical phonons and impurities and having 
in mind a systematic experimental thermal conductiv¬ 
ity study in the high-quality, single-crystalline, large- 
J spin-chain cuprates Sr 2 Cu 03 and SrCu02 that has 
been recently conducted PTHSO] . Several attempts to de¬ 
velop a suitable formalism to describe this phenomenon 
have been made in the past However, these ap¬ 

proaches either relied on unrealistic choices of parame¬ 
ters [2ll [26] or offered only qualitative insights [2ll ES] ■ 

Below we attempt to bridge the gap between experi¬ 
ment and theory. We argue that the heat conductivity by 
spin excitations can be quantitatively described within 
the bosonization framework with the large-momentum 
scattering by optical phonons or impurities. For weak 
impurities, scattering grows stronger at lower temper¬ 
ature, a feature intimately related to a critical char¬ 
acter of the 5 = 1/2 Heisenberg chains [21]. Taking 
into account multi-spin-boson processes, it follows nat¬ 
urally from our microscopic calculations that scattering 
by phonons bears a close similarity to that by weak impu¬ 


rities, except that the phonons are thermally populated 
and thus control heat transport at high T. This is also 
in accord with a physical picture of phonons playing the 
role of impurities for the fast spin excitations. Within 
this picture, the transport relaxation time is the same as 
spin-boson scattering time and the corresponding mean- 
free path fits excellently the available experimental data. 
Further systematic extensions of our theory to include 
multi-phonon scattering that can influence thermal con¬ 
ductivity at higher temperature are briefly discussed. 

Finally, we emphasize an important physical constraint 
on the strength of spin-phonon coupling of magnetoe¬ 
lastic nature EH isa, which is weak in the materials 
of interest. While an estimate of this coupling can be 
made microscopically, a simple piece of phenomenologi¬ 
cal evidence for this criterion is the absence of the spin- 
Peierls transition in real compounds down to very low 
temperatures. Our theory easily satisfies the proposed 
constraint, setting itself apart from the previous ap¬ 
proaches |2H US]' We thus provide a microscopic, in¬ 
ternally consistent description of thermal transport and 
scattering in ID spin chains, which satisfies naturalness 
criteria by having weak spin-phonon coupling and con¬ 
forming to an analogy between phonon and impurity 
scatterings. 

Spin-phonon coupling Hamiltonian .—The nearest- 
neighbor Hamiltonian of an 5 = 1/2 Heisenberg chain 
magnetoelastically coupled to phonons is 

(U) 

where (ij) denotes nearest-neighbor lattice sites. A stan¬ 
dard Jordan-Wigner transformation with the subsequent 
bosonization m and the lowest-order expansion in lat¬ 
tice displacements brings it to the following form: 

H = ^kh\bf. -\- "Hs-ph ) (2) 
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FIG. 1: (Color online) Schematics of the spectra of bosonic 
excitations in a large-J, 5= 1/2 Heisenberg spin chain (dis¬ 
persive branch Ek = nlkl and a continuum at Q = tt) and of 
the dispersionless optical phonon branch ljq (the horizontal 
line). 

where 5^,^^ represents spin-boson operators of the ex¬ 
citation with Efc = v\k\ (sketched in Fig. [^, velocity 
is V = 7rJa/2, k is the ID momentum, and a is the 
lattice spacing. Hamiltonian 'Hs-ph describes a large- 
momentum, q Ks Q = tt/o, spin-boson scattering by 
phonons 

'Hs-ph = J dxlJj;{x)cos(^^{x) + Qx^ , (3) 

where A = adJ/dx, x is the direction along the chains, 
the lattice displacement field 'U(x) is associated with the 
optical and zone boundary phonons, and the spin-boson 
field ^(x) = e^^^(bl + b_y.)/L\k\, in which L is 

the linear size of the chain and we used the Luttinger- 
liquid parameter /C = 1/2 for the Heisenberg case [33]. 
Small-momentum scattering is deliberately ignored, as 
the corresponding vertex carries small in-plane momen¬ 
tum of the phonon and leads to negligible scattering ef¬ 
fects [26] . 

We note that boson-boson scattering cannot dissipate 
the heat current [HKiiiiss] and thus is neglected. 

Self-energy and relaxation rate. —Assuming the spin- 
phonon coupling to be small, a conjecture discussed be¬ 
low in detail |32j , one can consider only the second-order 
spin-boson self-energy in Fig. [^a) given by 
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FIG. 2: (Color online) Multi-boson diagrams contributing 
to the scattering rate of spin bosons on (a) phonons and (b) 
weak impurities. Shaded ellipses represent a set of diagrams 
involving an arbitrary number of spin bosons in the interme¬ 
diate state. Solid and wavy lines are the Green’s functions of 
spin bosons and phonons, respectively. 


we obtain D{t,x) = a5{x)Fr{uio)/2muio with 


FriuJo) = noe‘^°'"{nol)e , (5) 

where np = 1/ (e“°/^ — 1) is the phonon distribution func¬ 
tion, m is the mass of the unit cell, and h=kB = i.. 

For the bosonic field correlator in the spin-phonon self¬ 
energy @ in Fig. I^a), we note an immediate similar¬ 
ity to the second-order T-matrix for the weak impurity 
scattering in Fig. [2|b), which also generates a large- 
momentum transfer [2^. The correlator can be evaluated 
at a; —0 and T <^J [M] |33] and leads to 


/g-i<I>(0.0)gi$(T,0)^ 


ttT 

J\ sin(7rrT)| 


( 6 ) 


Then the self-energy at Matsubara frequency is 


Alfe(w„) — 5gp 


2TJ e--- - 1 

a\k\ Jo I sin(7rrT)| ^ ’ 


(7) 


where we introduced a naturally appearing dimensionless 
spin-phonon coupling constant (?sp = X/{aJy/2muJo) [ST] 
122] . For the spin-boson scattering rate, we need the imag¬ 
inary part of the self-energy that is analytically continued 
to real frequencies. The transformations allowing us to 
perform the integration in Eq. 0 exactly are delegated 
to Ref. [33] ■ Here, we simply list the answer, 

ImEfc(a;) = ^ (2no -f l) (l - /+ - /_) , (8) 

where f± = l/(e‘^^‘^“ -I- 1). The fermionic distributions 
can be seen as a result of a re-fermionization of bosons 
via a multiple-boson scattering. The result ([^ can be 
expanded in w/T, yielding 


where D{t, x) = (U2,(0,0)U2,(r, x)) is the phonon propa¬ 
gator and the second quantization of lattice-displacement 
field is standard [36] . We exploit the large value of J com¬ 
pared to a typical Debye energy (in cuprates J/0 _d ~ 10), 
which allows us to neglect dispersion of the phonon 
branches near the 7r-point in Fig. [^ Then, the lattice- 
displacement correlator is fully local in space [33] and 
separates into a sum over phonon branches £ that have 
non-zero projections of their polarizations, on the 
chain axis x. Considering for simplicity only one lon¬ 
gitudinal phonon with the energy ujq (see Fig. [^, and 
reserving the right to add more phonon branches later. 


ImEfe(a;) 


2 2Ju} 1 

a\k\T sinh(wo/r)’ 


(9) 


which holds exceptionally well for all w < T of interest. 
Generally, the single-particle scattering rate ^ should 
differ from the transport relaxation rate, but for the 
impurity-like scattering the two become the same. 

Mean-free path. —Then, the on-shell approximation, 
w = £k, in Eq. 0 yields the inverse spin-boson mean- 
free path, Xjl=\lvT., due to spin-phonon scattering: 


/^A \ 2 1 

\ a / T sinh(wo/T) 


( 10 ) 
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This result is fc-independent and thus can be compared 
directly to the transport mean-free path extracted from 
thermal conductivity data m [29]. We note that the 
l/T prefactor in Eq. (10 1 is strongly reminiscent of the 
result for the scattering on weak impurities US] [371 |3B| : 
(£imp/a) ^ =»T'imp {•J/'T), where riimp is the con¬ 

centration of such impurities and SJ is the strength of 
impurity potential. Clearly, this scattering gets stronger 
with lowering T, down to the Kane-Fisher scale, Tkf oc 
SJ^/J, below which weak impurity becomes a strong 
scatterer, similar to a chain break |33|. This behavior is a 
consequence of a critical character of spin chains mEoi. 
Since phonons should be seen as weak impurities by the 
fast spin excitations, it is natural that the spin-phonon 
scattering yields the same l/T prefactor in Eq. (10). 

While the other thermal factor in (10), 1/sinh(a;o/T), 
does not coincide with the phonon population rig, both 
have the same high- and low-T asymptotes. For T <C 
wg, the mean-free path ( 10 ) exhibits activated behavior, 
similar to the findings of other works [24l I25j . 

In addition to the considered scattering mechanisms, 
the low-T spin thermal conductivity in real materials is 
limited by strong defects that act like chain breaks m 
I2H1I3S|. The corresponding mean-free path is an average 
length of a defect-free chain segment, l/£b = where 
nb is the concentration of these defects. 

Comparison with experiments. — Figure[^shows the T- 
dependence of the mean-free path of ID spin excitations 
in Sr 2 Cu 03 and SrCu 02 dTHMj- The data are extracted 
from the thermal conductivity measurements via a ki¬ 
netic relation, £{T) = K{T)/vCviT), using theoretical val¬ 
ues [41] for the specific heat of spin chains Cv{T) [oc T 
at T ^ J]. Because of high purity, the mean-free path 
exceeds lO^a at low T, with the difference between the 
two compounds due to the residual concentrations of the 
defects. The two sets of data become quantitatively very 
close at higher T, implying that a similar scattering is 
dominating propagation of heat in both materials [29j . 

Figure]^ shows our successful fits of the data by com¬ 
bining spin-phonon ( 10 ) and strong-impurity scatterings, 
via Matthiessen’s rule We note that 
the low-T part of the data, T < 40 K, has a large uncer¬ 
tainty due to the subtraction of the phonon part of the 
thermal conductivity (see [57| |in])j and that it can be 
fit with an equal success by a combination of weak and 
strong impurities, ^ ~-^bi^p+^b ^ • Since it is a secondary 
issue for our study, the simplest account of this regime 
by strong impurities suffices. To fit the spin-boson mean- 
free path above 40 K, we assume that the spin-bosons are 
scattered by two phonon modes with wgq = 300 K and 
<^ 0,2 = 740 K. Of the two, the first roughly corresponds to 
the longitudinal zone-boundary phonon and the second 
to the high-energy stretching mode |42l [43] , both likely 


having the strongest coupling to spin chains. In Eq. (10) 
we used the value of J = 2600 K [33] and the spin-phonon 
coupling constants gi^sp = 0.020(1) and 32 ,sp = 0.10(1) 



FIG. 3: (Color online) Mean-free path of spin excitations in 
Sr 2 Cu 03 and SrCu02 [371133] (the symbols). Lines are theory 
fits, see the text. Concentrations of strong impnrities, ng, 
phonon energies, wo.i, and spin-phonon coupling constants, 
indicated in the graph. 


have provided the fit in Fig. [^for Sr 2 Cu 03 and SrCu 02 . 
By choosing different wg^^’s, one can obtain somewhat 
different values of the 5 i,sp’s needed for the fit, but they 
never exceed or even reach the physical bounds discussed 
next. 

Bounds on spin-phonon coupling .—We now discuss 
physical bounds on the spin-phonon coupling constant 
<7sp = A/(JaA/2mwg). As discussed in Refs. lEI^, 
the constant is a product of two parameters, with one 
characterizing the change of J by atomic displacement, 
7 = X/J = a{dJ/dx)/J, and the other representing an 
amplitude of zero-point atomic motion relative to a lat¬ 
tice constant [36] . a = where h is made 

explicit and m is the reduced mass associated with the 
phonon mode wg. Parameter a is small, while 7 can be 
large [SI [33] because the superexchange is sensitive to 
the interatomic distance. Regarding cuprates, one can 
estimate a « 0.01. The superexchange parameter has a 
larger uncertainty, with indirect studies giving a range of 
7 = 3—14 [33] [ 35 ] and a consideration of a wider class of 
materials suggesting an upper limit of 7 <20 [S]- Thus, 
the microscopic upper bound on the spin-phonon cou¬ 
pling constant in ID cuprates can be put at g^‘^^^0.2, 
justifying the weak-coupling treatment of the spin-boson 
scattering on phonons in Eq. ([^. 

A less restrictive, but purely phenomenological crite¬ 
rion limiting the strength of the spin-phonon coupling is 
the absence of the spin-Peierls transition in ID cuprates 
down to about 5 K («0.002J), where the 3D Neel order¬ 
ing can be argued to preempt other transitions. Using 
Tp « this can be translated to the upper limit 

on the spin-phonon coupling 5 ™^’'« 0.35. 

We now offer a critique of the previous considerations 
of thermal transport in ID spin chains. In particular, in 
experimental works [i7|[ini[3i, the spin-phonon mean- 
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free path is repeatedly fit by the form = ATe~^ 

with w* r:! 200 K, inspired by the phonon-mediated Umk- 
lapp scenario [2l]- First, most of the data in Fig. [^should 
be outside the quantitative accuracy range of this ex¬ 
pression, which is limited to T<a;*/3«70K, as the ex¬ 
ponent is only a low-T limit of the phonon distribution 
function. More importantly, translating the values of A 
used in Refs. mM m] to the dimensionless coupling 
constant via A = g^^/Ja gives g^-p « 1, which is exceed¬ 
ingly large for the perturbative treatment to hold and 
lies way outside the allowed range. This strong coupling 
also implies a spin-Peierls transition at Tp ~ J, while 
no such transition is observed. Likewise, our previous 
study, which considered small-momentum scattering on 
acoustic phonon branches |26j . required an anomalously 
strong spin-phonon interaction gsp > 1. Thus, there is a 
serious “naturalness” problem with previous theoretical 
considerations. 



FIG. 4: (Color online) Same as in Fig. The solid line 
includes T-dependence in the spin-phonon coupling; see the 
text. 


In the present work, the dimensionless spin-phonon 
coupling constants are well within the range of the mi¬ 
croscopic expectations for the cuprates, (?sp = 0.02 —0.1, 
implying extremely low spin-Peierls transition tempera¬ 
tures. While the offered analysis of the physical bounds 
is not a proof of the validity of our theory, it is certainly a 
strong argument against the validity of the previous ap¬ 
proaches, which require unphysically large spin-phonon 
coupling. 

Multi-phonon scattering .—We note that for T > coq 
the spin-boson mean-free path in (101 saturates at 
(4p/a)~ «gsp2J/a;o- While this is not unphysical, one 
can still expect that the other, T-dependent terms may 
become important for T > ujq. Corrections of the or¬ 
der T/J are neglected in our derivation (see [33]), since 
Tj J is small in the relevant temperature range. Another 
possible source of the T-dependence is the multi-phonon 
scattering. Superficially, the two-phonon scattering pro¬ 
cesses have to be negligible because of the smallness of the 
spin-phonon coupling discussed above. However, there 
are factors that can compensate for this smallness. First, 
the two-phonon scattering is less restrictive, as the trans¬ 
verse phonons can also contribute. Second, in the non- 
Bravais lattices, the two-phonon processes are also am¬ 
plified by the number of atoms in a unit cell, Na- That 
is, for the single-phonon processes, the number of lon¬ 
gitudinal phonons that couple to spins via ([^ is Alj, of 
which we have chosen only two for our fits in Fig. On 
the other hand, when a spin-boson scattering is due to 
the emission or absorption of two phonons, the number 
of possible processes can be as large as 0{N^). A naive 
and certainly overly optimistic estimate of their number 
assuming independent polarization and a branch index 
for each phonon involved in the scattering yields (3A"a)^- 
In cuprates [331 j the total number of phonon modes is 
large, so this combinatorial factor can be substantial. 

A somewhat tedious, but straightforward algebra [33] 


yields the following result for the two-phonon scattering 

^ „ ./ rriQn ( / / l 

( 11 ) 


-sp,2 


2 J cosh(a;o/T) 

lisp,2 


T sinh^ (wq/T) ’ 


where gt^ oc C' 25 sp. When compared to (101, the re¬ 


sult in (11) contains an extra factor, gl ^04)1, and a 


large combinatorial factor, C 2 . Clearly, at T <^ujq, the 
two-phonon mean-free path follows the same behavior as 


(10), thus simply renormalizing single-phonon scattering. 


However, at T > wq, it carries an extra power of T/wq, 
( 4 p, 2 /a) ^ « ^sp ^d^T/wg, thus amounting to an expan¬ 
sion in TjujQ^i, which can be argued to follow naturally 
from the multi-phonon scattering processes. 

Without going into non-generic microscopic consider¬ 
ations, one can suggest a simple ansatz to account for 
the T/wo_i-expansion with the T-dependence of the spin- 
phonon coupling in the form ( 7 sp,i(T) = ^sp.i (l+?’i’T-o,i)j 
where no^i = l/(e‘^“''/’^ — 1) as before. This form meets 
both the low-T and the high-T behavior of the two- 
phonon mean-free path discussed above. A fit of the 
SrCu 02 data using this ansatz with = 1 is provided 
in Fig. The bare spin-phonon coupling constants gi^sp 
are even smaller than in Fig.[^ especially for the higher- 
energy mode. The result with the bare gasp’s is pro¬ 
vided for comparison. Although this figure is an illus¬ 
tration showing that our theory allows for systematic 
extensions by including multi-phonon processes, it also 
demonstrates a potential role of the latter in the T > wg 
regime and thus contributes to the general description of 
the heat transport in spin-chain materials. 

Conclusions .—We have provided a consistent micro¬ 
scopic theory for thermal transport and scattering in ID 
spin chains which stands out from previous attempts at 
such a theory by having weak spin-phonon coupling and 
conforming to the analogy of the phonon scattering to 
that on impurities. We have successfully fit the available 
experimental data and discussed possible extensions of 
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our theory for higher-T. Our approach should be appli¬ 
cable to the thermal conductivity in spin-ladder materi¬ 
als and can be extended to the transport phenomena in 
a variety of Luttinger liquids and ultracold atomic gases. 
Numerical verification of our results is also called for. 
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SPIN-BOSON SCATTERING ON PHONONS 


Phonon correlation function 


Bosonization and coupling to phonons 

Bosonization of the Jordan-Wigner fermion operators 
for the S' = 1/2 chain proceeds via introduction of the 
bosonic fields $ and 0 for the right- and left-moving 
fermion operators [T] according to 


Using standard quantization of the displacement field, 


phonon correlation function in (S5) is 


(U,(0,0)U,(t,R)) = -^ 


,-zPR 




N ^ 2mu;ot 


(S6) 


V'R = 


'0L = 


1 




na 


•\/27ra 


iv^(-3>-|-e) 


(51) 

(52) 


X ^^ap^(0)ap^(T)^ -|- ^a_p^(0)alp^(r) 


Assuming, for simplicity, that the dynamical force matrix 
at a given momenta can be diagonalized so that one of the 
polarizations (longitudinal) is along the chain and two 
remaining ones are perpendicular, then all of the relevant 


The bilinear combination of the fermionic operators polarizations yield (^p^) — 1 and the summation in (S6) 
V'rV'l be expressed as follows 




— i-\/27T^ 


2'Ka 


(S3) 


where the rescaled bosonic field $ = (1/-\/A)$ = 

The Tomonaga-Luttinger parameter JC = 1/2 corre¬ 
sponds to the Heisenberg limit of the model. 

Coupling to phonon field that leads to the large- 
momentum scattering can be written as 

Hs-ph = A / + h.c., (S4) 

where A = adJ/dx and \5x{x) is the a;-component of 
the lattice displacement field U for the optical and zone 
boundary phonons. 


Spin-boson self-energy 

Since we are interested in the scattering involving 
momenta close to the ID zone boundary momentum 
Q = 7r/a, one can assume that the corresponding phonon 
has zero velocity. Then we can write the spin-boson self¬ 
energy as 


Efc(r) = - 


2X^ 


■ J dxe^’^^ {Ux{0,0)Vx{t,x)) (S5) 

X \ c fc- 


will be over Na independent longitudinal phonon modes 
with individual frequencies where Na is the number 
of atoms per unit cell. 

If one were to include dispersion of phonons into con¬ 
sideration, the phonon propagator in (S6) will not be 


strictly local. However, for the fast spin excitations, the 
phonon does not propagate far during the time of inter¬ 
action. Thus, the effect of the phonon dispersion is small 
by the small parameter Qd/J ~ 0.1 and our u}o,t cor¬ 
respond to the average energies of the phonon branches 
that couple strongly to spins. 

Consider coupling to an individual phonon mode with 
the energy uiq. For this case we derive 


(U,(0,0)U,(r,R)) = 

X 

with uq = 


VoSjR) 

2inu}Q 

X [noe“«^-f (no-b l)e-‘^«^] 

1 


(57) 

5 

(58) 


where Vq is the volume of the unit cell, no is the phonon 
occupation number, and we assumed that r > 0. Since 
we need to evaluate the phonon propagator on a single 
chain, one can write Vo'^(R-) which leads to 


(U,(0,0)U,(r,a;)) 

a6{x) 

2mujQ 


(S9) 


[noe^^o" + (no + l)e-““"] 


'iTa^\k\ 


Obviously, for the case of Na phonon modes, one needs 
to reinstate the summation over them. 












Spin-boson correlation function 

We now turn to the spin-boson part. We have previ¬ 
ously obtained the spin-boson correlation function in [2] . 
Here we provide a more rigorous derivation of it. First, 

-iv'2Fi>(0.0)giv^4(T.0)\ _ g-^([^(0,0)-i>(r,0)]^^ (SIO) 


(Sll) 


Thus, we are to calculate g{T) defined by 
g(r) = -/[$(0,0)-$(r,0)''' 


= 2 ^$(0, 0)<i>(r, 0) ) - 2 ( $(0,0)<I>(0, 0)) . 
Therefore, we have 




L 


k^O 


OJ^ -b ’ 


(S12) 


with ujn = 27rTn, n = integer. 


Discarding the vanishing w„ = 0 term in (S12), we have 
the following identity for any a;„ ^ 0 


2T 


9(-) = -^ E -1) E 


OJn^O 


k^O 


oj'i -b 


(S13) 


= 2r E /*“ 




27r oj'i -b 


Here the exponent e is used to set an ultraviolet 

cutoff. Next, we rewrite 5 (t) as 


, , T f+^ dy ^ . 

1^ ■ 

.7-00 


- 1 


(S14) 


where y is the dimensionless integration variable and we 
used V = 7rJa/2. The summation over a;„ can be per¬ 
formed with the help of an identity 


E 


^uJri{ir-ylJ) 


UJri 


UJn>0 

As a result we obtain 
^ ^ I f dy 

X In 


In (l - 




1 -b cos(27rTT) 


dy 


271^ J 1 -b ?/' 


■ In 


sy 


sy"^ -b 4sin^(7rTr) 


(S16) 


where s = {2'kT/ JY 1. The remaining integral can be 

found in [3], which gives 


gij) = - In 




TT 2| sin(7rTT)|-b-^/s 


(S17) 


Exponentiating ( 7 (t), we arrive to the following expres- 

(S18) 


sion for the spin-boson correlator in (S5) 




2t:T 


ttT 


2 J| sin(7rT'r)|-b 27rr J| sin(7rrr)| ’ 

where in the last expression we neglected contributions 
of order 0{T/J). 


Evaluation of the self-energy 

Using the phonon propagator (|^ and the spin-boson 


correlation function (S18l obtained above and transform¬ 


ing the spin-boson self-energy Efc(r) in ( [S^ to the Mat- 
a frequer 

Efc (cOn) 


subara frequency domain yields Ek(u!n) given by 


r/3 giw„T _ ^ 

^ Jo ^'^|sin(7rrT)| 


\2ma^uJo / \ Ja\k\J 

[noe‘"“" + (no + 


(S19) 


Therefore, we need to evaluate the following integral 

I'd pil^nT _ 1 

4(wn) = / dr 

Jq sm(7rTr) 

for e = iwn. Here we used that t > 0. 


(S20) 


Consider Eq. (S20l for e = wq. Defining the complex 
variable 2 = ’^6 transform into an integral 

along the unit circle in the z-plane. This contour, how¬ 
ever, is not closed, since the circle is cut at z = 1 by a 
branch-cut running along the real axis from zero to -boo. 
Completing the unit circle to the contour C, see Fig. 
one can verify that 


la, 


r = — 

ttT 


P 1—7" 

dz- -0-*^o-i/2 = o, (S21) 

Ic k — z 


where ipo = wo/27rT, and contains two segments of 
the contour C along the real axis. 

Thus, Iipg reduces to the integral of the real variable 




dx 


1-x” 
1 — X 


-iipo-l/2 


. (S22) 



FIG. 5: (Color online) Contour C for the evaluation of 
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Method I 


The integral in (S22) can be evaluated using m , giving 


luio — 


1 + e' 


1^0 IT 


ttT 


V' ( n — i(po + 


-i) -iipo + 


(S23) 


where 'ijj{'w) is the digamma function. 

Ultimately, we need the real-frequency propagator. To 
obtain it, let us perform the substitution —>■ w -|- iO, 

or, equivalently, n —)■ ui/i2'i:T 0. In the real-frequency 
domain the following holds 


pU^olT 


Iljo ~ 


ttT 


ijj 


UJ 


Wo , 1 


\i2TrT i2'KT 2 


-V' 


Wq 

i2'nT 


1 


(S24) 


Imaginary part of can be expressed in terms of ele¬ 
mentary functions, since 


Im'i/) 



TT 

2 


tanh(7rj/). 


(S25) 


The latter equality can be derived with the help of the 
reflection formula for the digamma function '!/)(! — x) — 
ip{x) = 7rcot(7ra:), in which one has to substitute x = 
1/2 — iy. Consequently, 


Im 


1 -t 


2T 


tanh 


a; ± Wo \ 
2T J 


=p tanh 



(S26) 


Method II 



FIG. 6: (Color online) 


Contour C for I^ig 


in (S28l. 


standard trick utilizing the contour in Fig. Both yield 
identical results, coinciding with Eq. ( |S26[ ). The quantity 
Im I-gjg can be obtained from the latter result by simple 
substitution wq —— wq. 


Scattering rate 


The imaginary part of the self-energy in (S19) is pro¬ 


portional to a combination of Im I±uig with the phonon 
distribution functions that correspond to the phonon ab¬ 
sorption or emission processes 


F(w,wo) = noIm /,^0 -|- (no -I- l)Im/_,, 


(S29) 


where ng is defined in Eq. (S8). The expression for 


E(w, Wo) can be rewritten as a combination of the phonon 
and effectively fermionic distribution functions 

F(w,wo) = ^(l + 2no) (!-/+-/_), (S30) 


with f+ = —i— 

fjUJdzujQ 


-hi 


(S31) 


All of that finally yields the imaginary part of the (re¬ 
tarded) self-energy 


A substitution x = e transforms Iuj„ in (S22) to 


IbJg 


(1 -h e‘^«/'^) 


ttT 


• „-2nz _ 1 

dz ——-, 

sinh z 


(S27) 


Since we are, ultimately, interested in the imaginary part 
of Ilfc(w), it is the imaginary part of Ia;„ that we are 
concerned about. This is also a point at which an ana¬ 
lytical continuation can be made via iw„ —)■ w -I- iO, or, 
equivalently, n —>■ —iLij/2'xT. Introducing the variable 
Lp = oj/2 'kT and extending the limit to —oo we get 


ImSr(w) = -5. 


2J 

a\k\ 


(l + 2no)(l-/+-/-),(S32) 


where the dimensionless spin-phonon coupling constant 
5 sp = A/(aJ-\/2TOWo) is introduced. 

We would like to point out that the only approxima¬ 
tion that has been made in the derivation of our answer 


in (S321 from the Matsubara self-energy in (S5) is the 


omission of the 0{T/ J) terms in Eq. (S18). 

At w < T one can approximate the imaginary part of 


the self-energy in (S32) as 




(1 -h 

' ttT ' 


(S28) 


X 


sin(<^z) cos((2(/?o + ^)z) 
- 


The resulting integral can be found in or evaluated by 
using integration in a complex plane with the help of a 


ImE/®*(w) « -y. 


2 Jw 


a\k\T sinh(wo/r) 


(S33) 


We have verified the validity of approximation in (S33) 


by a numerical calculation of the relaxation rate using 


(S32) directly. The results differ by a slowly varying T- 


dependent function with the values near unity. 
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Finally, the scattering rate of spin-bosons on phonons 
is obtained by taking ImE™*(w) on the mass surface, 
w = v\k\, 


plph 

^ k 




TT,P 

~T~ 


1 

sinh {oJo/T) 


, (S34) 


where v = TTJa/2 was used as above. This single-phonon 
scattering rate vanishes exponentially as a function of 
temperature for T < ujq and saturates at a constant value 
for T > luq. 


TWO-PHONON SCATTERING 


Here we extend the formalism developed above to the 
interaction of spin-boson with two phonons. In the limit 
of purely Einstein phonons this procedure is straightfor¬ 
ward. We model the two-phonon coupling as 

^ E /+ h.c.) , (S35) 

ct—x,y,z 

rP‘ T 

where = a^ —. (S36) 

drl 

Needless to say, this Hamiltonian represents a greatly 
simplihed version of the general situation as the coupling 
constant is independent of both the branch index 
I and polarization of a phonon. The square of the 
phonon held, which enters this Hamiltonian, equals to 


= Z ^ + (S37) 

q,i,q' 

x(alp +a_ ,p). 


According to the Wick’s theorem, the corresponding 
propagator can be expressed as a product of two single¬ 
phonon propagators 


(U^(0,0)U|(x, r)) = 2 (U„(0,0)U„(a:, r))" <5„;3.(S38) 


The propagator itself is gi ven b y Eq. (S9). The singular¬ 
ity of the form [a5(a:)]^ in (|S38|) can be resolved through 


the usual prescription: [a(5(a;)]^ = aS(x). Consequently 


(U2(0,0)U|(a:,r)) 


(S39) 


aS{x) 

2m'^ujQ 


[noe““" + (no + l)e-‘"“l"<5„^. 


We note that for a multi-atomic unit cell the latter ex¬ 
pression has to be multiplied by N^. In other words, the 
multi-atomic unit cell enhances the contribution of the 
multi-phonon scattering. 

Applying the same technical approach used in the 
single-phonon calculation above to the two-phonon case 


is quite straightforward. The self-energy is 

v 2 ph, ^ (A2ph)2 3a nT 

2-11, — — 


(S40) 


^ 7ra®|fc| 2m?LOQ J 

X \nQl2bJo 4" 2no(^o + l)7o + (no + I)^/_2(jo] > 


where I± 2 luo,o are defined according to Eq. (S20). Then, 
one can obtain the scattering rate from this expression 

r,„, ^ co.h(.„/2r) 

Jm^a^ujQ T sinh^ (uJo/T) 


Using v = 'KaJ/2 and introducing dimensionless spin-two- 
phonon coupling constant, 2 = J)^/(2ma^a;o)^, 

we simplify the last expression to 


p2ph 


« g, 


J'^ cosh(a;o/2T) 
’ sinh2(wo/r) 


(S42) 


One can relate the two-phonon to the one-phonon cou¬ 
pling constant, ffsp ,2 C' 2 gsp, where C 2 is a large combi¬ 
natorial factor due to multiple phonon modes that can 
be involved in the scattering. At T <u!o, the two-phonon 
scattering rate yields the same exponential behavior as 


(S34). However, at T>uJo, it carries an extra power of 


T/wo 


p2ph 


2 J T 

^ 9sp,2 2 ’ 

Lc'n 


(S43) 


thus giving a natural expansion in T /wq in the multi¬ 
phonon scattering of spin-bosons. 


COMMENT ON PRIOR WORK 


The spin-phonon scattering considered in Ref. [21 while 


similar to (S4), included coupling only to acoustic 


phonons, which, in turn, necessarily lead only to a small- 
momentum scattering of spin excitations, with the cor¬ 
responding scattering rate carrying a smallness of a high 
power of T/J ^ 1 compared with the large-momentum 
scattering of the present work. As a result, the mecha¬ 
nism of [2] required large spin-phonon coupling constant, 
<7ac ~ 10, to explain experimental data. If the realistic 
coupling is used, the scattering processes studied in [2] 
turn out to be much weaker than those of the present 
work for most temperatures. 

However, at low temperatures, T <?; Odi the mecha¬ 
nism discussed in [2] can still be relevant as the expo¬ 
nentially small population of optical phonons in (S34) 


competes with the power-law smallness of the small- 
momentum transfer mechanism of [2]. In the notations 
of the present work, the scattering rate due to acoustic 
phonons at low-T can be written as [2] 


oc 5a 


2^5 

J40 


D 


(S44) 
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where g^c is the dimensionless coupling constant as de¬ 
fined in the present work. Making assumption on cou¬ 
pling constants being the same and Qd = <^o, one can 
compare this result with (S34) to find the temperature 
T* at which contributions of the two mechanisms to the 
scattering are equal. For the large-J systems, taking 
Qd/J « 0.1 leads to T* « O.10d. This estimate puts 
the regime at which the mechanism of [2] will be domi¬ 
nating to T < 30K, i.e. deep inside the regime which is 
controlled by impurity scattering even in the very clean 
materials discussed here. Thus, although, hypothetically, 
there is a regime where the mechanism of [2] will be dom¬ 
inating, it is mostly unimportant in the case of cuprates. 
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